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ABSTRACT 

Early -type galaxies (ETGs) define a narrow strip in the size -mass plane because 
of the observed correlation between the effective radius Reff and the total stellar 
mass Mi,. When expressed in logarithmic units, a linear relation, log R^ff oc 7logM^, 
is indeed observationally found, but the slope 7 deviates from the canonical 7 = 1/2 
value which can be naively predicted for a spherically symmetric isotropic system. We 
propose here that a transfer of angular momentum to the stellar component induces an 
anisotropy in the velocity space thus leading to a modified distribution function (DF). 
Assuming an Osipkov- Merritt like anisotropic DF, we derive an analytic relation be- 
tween the slope 7 of the size - mass relation and the slope a of the angular momentum 
term in the DF. With this simple model, we are then able to recover the observed 7 
value provided a is suitably set. It turns out that an anisotropy profile which is tan- 
gential inside 0.6ra and radial outside, with Va the anisotropy radius, is able to re- 
produce the observed size - mass relation observed for massive (M^ > 2 x 10^" ^q) 
elliptical galaxies. 
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1 INTRODUCTION 

Notwithstanding the quite large range spanned by their mor- 
phological and photometric properties, early -type galax- 
ies (hereafter, ETGs) show several interesting correlations 
among their colors, luminosities, velocity dispersions, effec- 
tive radii and surface brightness (e.g. Baum 1990, Faber & 
Jackson 1976, Kormendy 1977, Djorgovski & Davis 1987, 
Dressier et al. 1987, Bernardi et al. 2003a) . As a most famous 
example, ETGs are known to populate a tight plane, referred 
to as the fundamental plane, in the logarithmic space defined 
by the their effective radius Reff, intensity /e = I{Reff), 
and central velocity dispersion ctq (Djorgovski & Davis 1987, 
Dressier et al. 1987, Bender et al. 1992, Burstein et al. 1997). 
This fundamental plane is usually parameterized as : 

Keff oc (To le 

with (a, b) predicted to be (2, 1) if ETGs are in virial equi- 
librium. The observed plane is, however, tilted with respect 
to the virial one since the different determinations of {a,b), 
depending on the photometric band and the sample used, 
are always different from the virial values. Jorgensen et al 
(1996) first derived a = 1.24 ± 0.07 and b = -0.82 ± 0.02 



from a set of 225 early - type galaxies in nearby clusters ob- 
served in the r-band. While this result is consistent with 
the original observations of Djorgovsky & Davis (1987) and 
Dressier et al. (1987), it is nevertheless in striking contrast 
with the most recent determination relying on ~ 9000 ETGs 
observed within the framework of the Sloan Digital Sky Sur- 
vey (SDSS). Using this large sample, Bernardi et al. (2003a) 
have found a = 1.49 ± 0.05 and 6 = -0.75 ± 0.01, which are 
more similar to the K-band fundamental plane of Pahre 
et al. (1998). Although part of this discrepancy may be al- 
leviated taking care of the selection effects and the fitting 
method adopted, the reason for such an inconsistency is not 
currently clear. 

While the precise values of the FP coefficients are still 
debated, it is nevertheless clear that the observed FP is tilted 
with respect to the virial prediction. Such a tilt could be 
caused by a variation in the dynamical mass - to - light ra- 
tio for ETGs as a result of a varying dark matter fraction 
(e.g., Padmanabhan et al. 2004; Boylan - Kolchin et al. 2005) 
or stellar population variations (e.g., Gerhard et al. 2001). 
It is also worth noting that non - homology in the surface 
brightness profiles of elliptical galaxies (e.g., Graham & Col- 
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less 1997; Trujillo et al 2004) may be an other explanation 
of the fundamental plane tilt. 

Projections of the fundamental plane are also of interest 
in studies of galaxy evolution. The SDSS team (Bernardi 
et al. 2003b) measured both the radius - luminosity and the 
Faber - Jackson (1976) relations obtaining respectively: 

o r0.630±0.026 
r0.250±0.012 

Converting stellar luminosity in stellar masses according to 
the prescription of Kauffmann et al. (2003), Shen et al. 
(2003) then determined a size -mass relation as: 

for ETGs with M* > 2 x 10^° Mq. Finally, the Kor- 
mendy relation (1985) between effective surface brightness 
fj,e and effective radius has been measured as : 

He = (2.84 ± 0.06) log Reff + (15.647 ± 0.032) . 

It is worth noting that, although the FP itself is quite tight 
over a wide range of ETGs (e.g., Bernardi et al 2003a), there 
are a number of indications of variations in the FP projec- 
tions themselves with galaxy luminosity. Both the FP and 
its projections (and their eventual correlation with luminos- 
ity) provide strong constraints to any theory of formation 
and evolution of these galaxies. Scaling relations arc indeed 
often considered as major constraints for models of galaxy 
formation and evolution. For instance, they can be used to 
trace the contribution of the different stellar components to 
the total luminosity of galaxies, as in the case of the color - 
magnitude relation (sec, e.g., Visvanathan & Sandage 1977; 
Bower et al. 1992 for ETGs, and Tully et al. 1982, Gavazzi et 
al. 1996 for late - type galaxies) or to study the relationship 
between kinematical, structural and stellar population prop- 
erties of galaxies as in the case of the Tully - Fisher relation 
for spirals (Tully & Fisher 1977) and the fundamental plane 
for ellipticals. Moreover, the study of these different scaling 
relations has been crucial for showing the role of mass in 
the formation of galaxies (Gavazzi et al. 1996; Boselli et al. 
2001), a result now generally called downsizing effect, which 
is a new major constraint for hierarchical models of galcixy 
evolution (De Lucia et al. 2006) . 

A particularly interesting role has been played by the 
Kormendy relation in elucidating important differences be- 
tween giant and dwarf ellipticals. Dwarf ellipticals (dEs) are 
considered to be dark matter dominated systems and are 
therefore of vital interest in understanding galaxy formation 
in general (sec, e.g., Ferguson & Binggeli 1994 for a com- 
prehensive review). Low luminosity elliptical galaxies are 
distinguished from late -type galaxies by their smooth sur- 
face brightness profiles. For systems with luminosity smaller 
than Mb ^ —18, the smooth profile galaxies divide into two 
classes : compact galaxies with high central surface bright- 
ness (e.g. M32), and diffuse galaxies with low central surface 
brightness (e.g., the Local Group dwarf spheroidals). In the 
hierachical galaxy formation scenario dEs are supposed to be 
formed from average amplitude density fluctuations (Dekel 
& Silk 1986) where supernova driven galactic winds expel 
the gas content and reshape the galaxy body. In an alter- 
native view, dEs could form from progenitors through the 
process of galaxy harassment (Moore et al. 1998). Likely 



candidates for the progenitors could bo stripped late type 
galaxies. Data from the Local Group dEs and few galax- 
ies in the Virgo cluster (Bender & Nieto, 1990; Held et al. 
1990) indicated that dEs are supported by anisotropy thus 
introducing a dichotomy in the otherwise linear sequence of 
increasing rotational support with decreasing luminosity for 
the whole class of ellipticals. 

Furthermore, the study of the B band surface brightness 
vs. absolute magnitude relation and of the Kormendy rela- 
tion has originally shown a strong, apparent dichotomy in 
the behavior of dwarf and giant ellipticals. While in dwarfs 
the effective or central surface brightness increases with lu- 
minosity, the opposite trend is soon in giants (e.g. Ferguson 
& Binggeli 1994; Graham & Guzman 2003). An opposite 
trend between giants and dwarfs has been also observed in 
the Kormendy relation (e.g. Kormendy 1985; Capaccioli et 
al. 1992). This surprising result has been originally inter- 
preted as a clear indication that dEs are not the low lumi- 
nosity extension of giants, but rather an independent class 
of objects. The study of Graham & Guzman (2003), based 
on HST data, has however shown that this dichotomy is 
only apparent since it is due to a gradual steepening of the 
central radial profile with luminosity. On the other hand, 
several recent observational evidences and simulations seem 
to indicate that local group dwarf spheroidals (Mayer et al. 
2006), Virgo cluster dEs (Barazza et al. 2002; Conselice et al. 
2003; van Zee et al. 2004; Mastropietro et al. 2005; Lisker et 
al. 2006a,b; 2007, 2008, Lisker & Han 2008: Wichiclscn et al. 
2008) or generally dwarf spheroidals in other clusters such 
as Coma (Smith et al. 2008), Perseus (Penny & Conselice 
2008) or in the SDSS (Haines ct al. 2007) might be late- 
type galaxies recently perturbed by an hostile environment 
through ram pressure stripping or galaxy harassment. 

It is of peculiar interest the scaling relation connect- 
ing effective radius and mass in dEs and E galaxies. As no- 
ticed by Shen et al. (2003), this scaling relation in the case 
of E is Reff oc M°'^^, while in dwarfs the slope decreases 
to Reff oc M'^ '^*. The same behavior has been observed 
by Dabringhausen et al. (2008) thus arguing in favour of a 
different origin, such as the dE galaxies being tidal-dwarf 
galaxies (Okazaki & Tariiguchi 2000). If this were true, then 
dE galaxies and E galaxies, that presumably formed from 
a monolithic and rapid collapse, would have fundamentally 
different phase space distribution functions (DEs). One way 
to study a possible indication for a dynamically different 
origin of dE and E galaxies is to solve for a possible in- 
dication of different angular momentum dependency of the 
DF. With this contribution we perform a flrst step in this 
direction by assuming spherically symmetric models for dE 
and E galaxies, but allowing the DF to depend on the an- 
gular momentum. Angular momentum is acquired through 
tidal torques of the protogalaxy with the neighboring ones 
(Hoyle 1949; Peebles 1969; Hoffman 1986a,b; Eisenstein & 
Loeb 1995; Catelan & Theuns 1996a,b). Baryons and the 
halo particles are well mixed initially (Mo et al. 1998; Car- 
done & Sereno 2005; Treu & Koopmans 2002; Keeton 2001). 
This is due to the fact that the original angular momentum 
of the dark matter halo comes from gravitational (tidal) in- 
teractions with its environment. Thus, the dark matter and 
the gas experience the same torque in the process of halo as- 
sembly and should initially have (almost) the same specific 
angular momentum (Klypin et al. 2002). 
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In the following, we shall study how an injection of an- 
gular momentum in the system changes the scaling relation 
Reff -Mi,. In a spherically symmetric isotropic system, one 
expects that Reff oc MJ, with 7 = 1/2, while observations 
give 7 — 0.56 (Shcn ct al. 2003). Wc therefore investigate 
whether the injection of angular momentum may help in ex- 
plaining the observed deviation from the naive expectation. 

The paper is organized as follows. Section 2 describes 
the Sersic profile and the Prugniel and Simien model used 
to deproject it. Sections 3 and 4 are devoted to the calcula- 
tion of the model DF in the isotropic and anisotropic case, 
respectively. The resulting size -mass relation is discussed 
in Sect. 5, while we derive a similar relation for dark halo 
models in Sect. 6. Conclusions are then presented in Sect. 7. 



2 THE PS MODEL 

Notwithstanding the large range in luminosity, mass and 
size, ETGs represent a surprisingly regular class of objects 
concerning their photometric properties. Indeed, as well 
known (Caon et al. 1993; Graham & Golless 1997; Prugniel 
& Simien 1997), their surface brightness is well described by 
the Sersic (1968) profile: 



with 



I{R) = le exp • 



R 



R. 



1/n 



(1) 



with R the cylindrical radius* on the plane of the sky and 
le the luminosity intensity at the effective radius Ref / . The 
constant 6„ is determined by the condition that the luminos- 
ity within Ref / is half the total luminosity, i.e. hn is found 
by solving : 



r(2n,6„) = r(2n)/2 



(2) 



where T(a, z) is the incomplete F function. Although Eq.(2) 
is straightforward to solve numerically, a very good analyt- 
ical approximation is given by (Ciotti & Bertin 1999) : 

1 0.009876 

6„ = 2n - . 

3 n 

The deprojection of the intensity profile in Eq.(l) is 
straightforward under the hypothesis of spherical symme- 
try, but, unfortunately, the result turns out to be a some- 
what involved combinations of the unusual Meijer functions 
(Mazure & Capelato 2002). In order to not deal with these 
difficult to handle expression, we prefer to use the model pro- 
posed by Prugniel and Simien (1997, hereafter PS97) whose 
three dimensional luminosity density reads : 



3(r) = jo 



Reff 



exp 



R, 



eff 



1/n 



(3) 



* Note that we have implicitly assumed that the intensity I does 
not depend on the angular coordinates. Actually, the isophotes 
are not concentric circles, but rather ellipses with variable ellip- 
ticitics and position angles so that / = I{R, ip). However, in order 
to be consistent with our assumption of spherical symmetry of 
the three dimensional mass profile, we will neglect such an efl'ect 
and, following a common practice, circularize the intensity profile 
considering circular isophothes with radii equal to the geometric 
mean of the major and minor axes. 



loh. 



"(l-Pn) 



F(2n) 



^° ~2Reff F[n(3-p„)] • 

Here, Jo — I{R = 0) = lec''" , while the constant p„ is chosen 
so that the projection of Eq.(3) matches a Sersic profile with 
the same values of {n,Reff,Ie)- A useful fitting formula is 
given as (Marquez et al. 2001) : 



_ ^ Q 0.6097 ^ 0.00563 



n 

Because of the assumed spherical symmetry, the luminosity 
profile may be simply obtained as : 



L(r) = 47r 



/ r'^j{r')dr' 
Jo 



which, for the PS model, becomes : 

^ 7[n(3-pn),6n??^^") 

^("^-^"^ FK3-P.)] 

where the total luminosity Lt reads : 

Lt = 27rn6;;^"e'"'F(2n)/ei?e// ■ 

Note that the total luminosity is the same as the projected 
one for the corresponding Sersic profile as can be immedi- 
ately check computing : 



(5) 



(6) 



Lt = 27r 



f 

Jo 



I(R)RdR 



As a final remark, let us stress that, under the hypothesis of 
constant M/L ratio, we can convert the luminosity density 
j(r) in a mass density p(r) simply as p(r) = T*j(r) so that 
the total mass of the stellar component reads = T*Lt. 
This may be determined from the measurement of the pho- 
tometric parameters (n, Reff, le) provided that an estimate 
of the stellar Af/L ratio T* is available (for instance, from 
the relation between T* and the colors or from fitting the 
galaxy spectrum to stellar population synthesis models). It 
is, finally, worth stressing that, according to Eq.(6), we ex- 

1/2 

pcct a size -mass relation as Reff oc Af^, as a consequence 
of the total luminosity being the same £is the projected one. 
Actually, some deviations from this simple scaling may be 
come out if we give off the homology hypothesis, i.e. that 
the Sersic index n is the same for all the galaxies. Indeed, 
should n systematically varies with the total stellar mass 
Mi,, e.g. as n oc M^ , then Eq.(6) gives Reff oc Mi^~"^^^. 
However, the scaling n oc M^ is far to be verified so that we 
will assume that the homology hypothesis indeed holds, at 
least as a first approximation. 



3 THE ISOTROPIC DF 

Assuming isotropy in the velocity space, the distribution 

function (hereafter, DF) of a spherically symmetric model 
may be easily recovered resorting to the Eddington formula 
(Binney & Tremaine 1987) : 



fin = 



dp 



1 d r 

^■K'^dSj^ 

1 d^p d^ 



* ^/£\d'^!J^=o 



(7) 
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with £ — >]/ — 17^/2 the total energy per unit of mass and 
^{r) = — $(r) the gravitational potential (with the sign 
changed). Note that we will refer to the second expression 
since it avoids differentiating the integral which has often 
to be evaluated numerically. Moreover, the second term in 
this latter expression is typically negligible because the mass 
density scales faster than in most models. 

As a preliminary step, we have therefore to evaluate 
$(r) by solving the Poisson equation — AirGp. Due to 
the spherical assumption, this reduces to : 



1 d / 

dr V dr 



= 'iivGp 



It is convenient to introduce the following dimensionless 
quantities : 



ri = r/Reff , ^ = <i>/4:TvGpoR: 



with po = ^*jo so that the Poisson equation becomes : 



1 d 



Tj^ dr) \ drj ^ 



expl 



-6n(77^/"-l)] 



(8) 



Integrating first and imposing that the force F oc —d^/drj 
vanishes at the centre gives : 



d-q 



nb' 



-n(3-p„) 



7[n(3- 



(9) 



so that, after a second integration, we finally get : 



$(7,) = -n6-"(3-''"'7[n(3-p„),b„r, 
n6-"(2-''"'r[n(2-p„),b„r,i/" 



l/r 



(10) 



having imposed that the potential is null at infinity. We 
stress that the central scaled potential only depends on n : 

n6;;"<''-''"'r[n(2 - 



lim ${r]) 



(11) 



In order to use the Eddington formula, we should now invert 
the relation \l/ — >I'(r) and replace r(^) into Eq.(3) to get 
p('I') and then dp/d^ . Needless to say, this is not analytically 
possible so that we prefer to escape the problem by a simple 
change of variables. First, let us consider: 




(12) 



having trivially defined * = v|//47rGpo-Re// = — i, 5 — 
£ /'iirGpoR^f f , p = p/po, and used the property that 'if van- 
ishes at infinity. We can now proceed the same way to get : 




0.4 0.6 
efflO) 

Figure 1. The normalized isotropic DF f{£, n) as function of the 
normalized energy £/^{0) for n ranging from 1 to 9 in steps of 
2 (from the bottom yellow to the top red curve) with the black 
line referring to the de Vaucouleurs (n = 4) case. Note that, 
due to numerical reason, ^(0) is actually ^(10~*) without any 
significant loss of precision. 



(13) 



Inserting Eqs.(12) and (13) into the second of Eq.(7), the 
isotropic DF finally reads : 




fin 



G^/^pl^'^Rf 



(14) 



■<=// 



with the dimensionless DF defined as : 



+ 



29/2,^7/2 

dr; 

IN 3/2 1 




dp d^p 



/M- 1 rf£ ld_p 



drip J dr] drf 
^ ' dp 



with r) the (numerical) solution of : 
*(77) = £ . 



(15) 



(16) 



Using Eqs.(3) and (10), it is only a matter of algebra to com- 
pute the different terms entering Eq.(15) for the PS model. 
We do not report here the full expressions for sake of short- 
ness, but we stress that the second term identically vanishes, 
while the final DF is a function of the scaled energy £ and 
the slope n of the Sersic profile. Moreover, since putting 
together Eqs.(4) and (6) gives: 



Po = T*jo = 



,n{3-p„) 



47ri?3^.^ nr[n(3-p,0] ' 
we can rewrite the isotropic DF for the PS model as : 



^3/2^1/2^-1/2^3^^ 

with Jo = I{R = 0) = he- 



Hn)f{£,n) 
G^I'^mV^F^''' 



(17) 



(18) 



and : 
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Fig. 1 sliows the normalized DF f{£, n) as function of £ 
over the energy range *(r;min) < £ < "^{Vmax) with 
) = (10 10*). It is worth noting that the re- 
sulting DF is always positive that is to say the PS model is 
a physical one whatever the value of n is. 



4 THE ANISOTROPIC DF 

According to the model wc are investigating, the PS den- 
sity profile and DF describe the initial configuration of the 
baryons which will then originate the visible part of the 
galaxy. Let us now assume that a physical mechanism takes 
place injecting angular momentum into this initial configu- 
ration. Such an injection can be attributed, for instance, to 
the action of tidal torques between neighboring protogalax- 
ies. As the protogalaxy collapses, it spins more and more 
rapidly thus acquiring angular momentum. As a further pos- 
sibility, it is worth remembering that a merger of two disc 
galaxies also produces a spheroidal stellar system retaining 
angular momentum from the original rotating discs. Inves- 
tigating the details of how angular momentum is injected 
in the initial PS - like configuration is outside our aims here, 
but we just recall these two possible scenarios to show that 
such a transfer of angular momentum is indeed possible. As 
a consequence, the phase space distribution of stars will be 
altered thus deviating from the isotropic DF we have con- 
sidered above. 

For spherically anisotropic systems, the DF may depend 
not only on the energy £, but also on the total angular 
momentum L. In particular, we will focus our attention on 
the Osipkov-Merritt (OM) models (Osipkov 1979; Merritt 
1985) where the DF reads : 



f{£,L)=MQ)L' 



(20) 



where the isotropic DF /o is evaluated in the lowered bind- 
ing energy Q = £ — L^/2r^. The slope parameter a and 
the anisotropy radius Va determine the velocity dispersion 
anisotropy profile given by : 



-\- ri 



(21) 



For a positive a, the anisotropy is tangential for < ar^ 
to become then radial for r^/r^ > a. On the other hand, a 
negative a gives a model which is radially anisotropic every- 
where whatever the value of Va is. Note that, for a = and 
To = 1, we obtain a model which is isotropic in the inner 
region and radially anisotropic in the outer ones. 

The density profile corresponding to the anisotropic OM 
models may be computed as (Cudderford 1991) : 



p(r) = 



(27r)^/^(2r^)" V{a + l) 
(l + rVr2)«+i r(a + 3/2) 



(Q)(*-Q)"+^/'rfQ . 



(22) 



Inserting Eq.(18) for the isotropic DF and changing variable 
from Q to Q = Q/^nGpoReff, we easily get : 



P(*,f?) = 



(27r)^''-(2i^:^-^-)",/-" r(» + l) 
{1 + Vyvl)"+' r(a + 3/2) 



{4nGpoRljjr+'/' 



2(6a-H9)/2^a-H3p(^ -f 1 



(a+3/2) 



r(a + 3/2)*- 



w r'O: a+1 r>4a 



a+l 



f{q,n){^-qr+'/^dq (23) 



where we have also introduced the more manageable vari- 
ables V = */^o, q = Q/*o, with *o = *(?) = 0). 

As an important step, we have to check that the right 
hand side of Eq.(23) has indeed the physical dimensions of 
a mass density. It is straightforward to verify that this is 
not the case because of the term (GpoReff)"' ■ The origin of 
this discrepancy may be, however, easily understood noting 
that GpoRiff has the same physical dimensions as L^, i.e. 
the square modulus of the angular momentum. Actually, in 
Eq.(20), we have implicitly assumed that L is dimension- 
less in order to follow the common practice. However, more 
correctly, we must introduce an (up to now) arbitrary nor- 
malization constant angular momentum Lq to recover the 
correct result. Taking care of this and using Eq.(ll) for V&o, 
we therefore rewrite the mass density as : 



2(^'^+^)/^Tr"+^T{a + l) / nr[n(2-p„)] 



r{a + 3/2) 

Leff\'^ po»?'"p(^,n,a) 
+1 



, n(2-p„) 



V Lo J (1 + 7?V^2)" 
having naively defined^ : 

Leff 



»+3/2 



(24) 



(25) 



(26) 



p{ip,n,a)= f{q,n){ip - q)"^^^^dq 
Jo 

Eq.(24) gives p as function of xp and rj, while one is usually 
interested in p as function of rj only. To this aim, we have to 
first solve the Poisson equation for the scaled dimensionless 
potential tp. Starting from 

V^* = -47rGp , 

it is immediate to get : 

772 dri\ dri J ^opo ' 
so that we have to solve the following differential equation : 



v> = -- 



r(a + 3/2) 



,n(2-pn) 



\ Lo J (1 



(l + r;Vr?2)"+i ' ^^^^ 

Needless to say, solving this equation analytically is not pos- 
sible given its high nonlinearity because of the way tp enters 



t It is worth stressing tliat L^j j is not the value of the angu- 
lar momentum at the effective radius. Although not completely 
correct, an order of magnitude estimate of i/e// is given by 
■sjGM^R^f I, with Mt the total mass of the luminous component. 
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Figure 2. The shape of the density profile of the OM-PS model. In the left panel, we set log{L^f f /Lq) = 0.0 and a = 0.15 (short 
dashed), 0.50 (solid), 0.75 (long dashed), while in the right panel a = 0.50 and three values of \og,{L^f j / Lq) are considered, namely 
—1.0 (short dashed), 0.0 (solid), 1.0 (long dashed). We arbitrarily set n = 4 and rja = fa/Ri^ff = 0.01 in both panels. 



in the determination of p(V', n, a). However, a numerical so- 
lution is possible provided one sets the Sersic index n, the 
anisotropy parameters {a, rja), and the scaled characteris- 
tic angular momentum^ L^ff/Lo. Once such a solution for 
■)/; has been obtained, we can numerically invert the relation 
ij) = ipiri) to finally get p = p{ri). It is, however, worth stress- 
ing that, as a consequence of how ip{ri) has been obtained, 
the mass density pijj) will bo parameterized by the quantities 
(n,a,ria,Leff/Lo) which are therefore the four parameters 
needed to assign the anisotropic OM - PS model. 

It is interesting to look at the density profile of the OM - 
PS model. To this end, we show in Fig. 2 the normalized den- 
sity ps = Ti^^/il + V^/'na)"^^p{v) for models with n = 4 ar- 
bitrarily setting ria = 0.01. As a general result, we note that 
an inner core developes when adding angular momentum to 
the orbit distribution. As left panel shows, the higher is a, 
the more prominent is the core, i.e. ps ^ const for r] < rjc 
with r]c an increasing function of a. However, outside this 
inner region, the model resembles the original PS one so that 
it is likely that its projection still recovers the Sersic surface 
brightness profile with great precision everywhere but for 
1] < r]c. Since rjc is typically quite small {ric ~ 0.01), we can 
still use the value of n retrieved from fitting the Sersic law 
to the observed photometry, which is what we will assume 
hereafter. Concerning the role of Lef / /Lo, right panel shows 
that it is mainly a scaling parameter shifting up the full pro- 
file as Leff/Lo becomes smaller^. While changing n leaves 
qualitatively unaltered the above results, the value adopted 
for rja plays a most significant role. Indeed, depending on the 
{a,Leff/Lo), should rja be larger than a critical value, Ps{rj) 
can also become an increasing function of rj in the very inner 
regions thus leading to an unphysical density profile. Inves- 
tigating in detail what is the region of the 4D parameter 



■t Note that the actual value of Lq is up to now meaningless since 
what is needed to solve Eq. (27) is the ratio L^ff/Lg. Changing 
the value of Lq gives only a trivial rescaling of the results. 
§ Note that, in Fig. 2, wc have used log{Lgff /Lq) as variable 
so that we are actually investigating a range covering two order 
of magnitudes. For this reason, the inner core of the model with 
log {Lgff/Lo) = 1.0 shifts outside the plot, but it is still present. 



space {n, a, ria, Leff/Lo) giving rise to physically meaning- 
ful OM-PS models is, however, outside our aims here. As 
a final remark, we warn the reader to not use Eq.(17) to 
replace po with the total mass in Eq.(25) since Eq.(17) 
only holds when p{r) has the functional expression given by 
Eq.(3) which is not the case here. Nevertheless, one can still 
resort to Eq. (4) since this is an outcome of the requirement 
that the projected PS model matches the Sersic profile which 
is still approximately true for the OM - PS case (except in 
the very inner regions). 

5 THE SIZE -MASS RELATION 

Notwithstanding the anisotropy in the velocity space, the 
spherical symmetry of the system makes it possible to com- 
pute the total mass as : 

Mi, = 4n / r^p{r)dr = AirRlff / r]^ p{i])dri . 
Jo Jo 

It is now only a matter of algebra to insert Eq.(23) into the 
above relation to get : 

^ 2<^"+^^)/V"+"r(Q + 1) f e''"b;;'^-^"'r(2n) ) 
r(a + 3/2) \ r[n(3-p„)] J 

f nr[n{2-p„)] Y+'^' fLeffY\ r o2 

X I{n,a,r]a, Leff/Lo) , (28) 

where we have used po = T*jo with jo from Eq.(4) and 
defined : 

lin, a, Va, Leff/Lo) = ^l^^^+. drj . (29) 

Note that, when a = 0, the density profile reduces to the 
original PS one and we get the usual result M* cx Reff, i-e. 
logi?e// oc (1/2) log M*, as expected. 

Let us now suppose that, for a fixed value of n, the mass 
integral (29) can be approximated as : 

logT = a + ba + c\ogjCeff + d\og'na 
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+ ea log Ceff + fa log rfa 



(30) 



(a, b, c, d, e, f) ~ (-2.357, 1.780, -0.017, 0.223, -2.316, 0.114) 



with Ceff = Leff/Lo and {a,b,c,d,e,f) constant param- 
eters to be determined. Inserting this ansatz into Eq.(28) 
gives : 



M* = 



-tt" 'r(a - 1) j o' '"l);. ""'r(2n) 
r[n(3-p„)] 



nr(2n 



r(a + 3/2) 

a + 3/2 



X ^(a) , 

with : 

yt(a) = 10"+''" . 

Using Eqs.(4) and (25), it is trivial to show that : 



(31) 



^c+{e+2)a 



r(2n) 



2r[n(3-p„)] 

c+(e + 2)a 



c+(e + 2)Q 



-R. 



3[c+(e+2)al 
2 

ff 



(32) 



Inserting this relation into Eq.(31) and using logarithmic 
units, it is then only a matter of algebra to finally get : 



log Reff = log RlO + 7 log 



2x1010 /i-i M© 



with : 



(3e + 6)a + (3c + 4) 



log Rio = 7 log (2x lO^^/i^^Mo) - log A/'(a) 



(33) 



(34) 



(35) 



... . _ 2(«°+»)/ V"+^r(a + 1) [ nr[n(2 - p„) ^ ^ "^^^^ 

7V(Q;j — , Q/o^ 1 ,n(2-p„) 



r(a + 3/2) 5; 

(e + 2)a + (c+2) 

r[n(3-p„)] 

(e+2)c + c 



T.7e^(a) 



2Ll 



n. 



Ja+d 



(36) 



Should all the ETGs have the same values of the photo- 
metric parameters (n, /e), the stellar M/L ratio T* and the 
anisotropy parameters (a,r;a), then their sizes and masses 
should turn out to be perfectly correlated as predicted by 
Eq.(33). Observationally, such a correlation is indeed found, 
but with an intrinsic scatter which may be obviously due to 
the scatter in the quantities entering the definitions of 7 and 
RiQ. Assuming, however, as a first approximation that n = 4 
for all galaxies (i.e., all the ETGs follow the de Vaucouleurs 
surface brightness profile) and that (a, 77a) do not change on 
a case -by -case basis, we can then estimate a and Lq from 
the observed values of 7 and iiio once the fitting coefficients 
(a, b, c, d, c, /) have been determined. This latter task can be 
eEisily done by computing the mass integral over a grid in 
the {a,Ceff-,r)a) space for a fixed value of n. In particular, 
for n = 4, we find : 



fits very well the mass integral in Eq.(29) with a mean per- 
centage residual (AX/T) ~ 0.1% and {A.X/X)rms — 6%. 
With these values, the slope of the size -mass relation turns 
out to be : 

7 ~ (1.974 - 0.474 X a)"^ 

so that, setting 7 = 0.56 as observationally found by Shen et 
al. (2003) for ETGs with M* > 2x10^° /i"^ M®, we finally 
get a ~ 0.40, i.e. the velocity anisotropy is tangential for 
r/ra < \fa. — 0.63 and radial elsewhere. For lower mass 
systems, Shen et al. (2003) find 7 = 0.14 thus giving a = 
—10.2. Such a large negative value should argue in favor of a 
fully radial anisotropy, but also imply a quite strong angular 
momentum term in the anisotropic DF which seems quite 
unusual. Moreover, we have computed the mass integral only 
for Q > so that our approximation for log! should not be 
extrapolated to models with negative a. 

It is worth stressing, however, that the surface bright- 
ness profile of dwarf ellipticals is better fitted by Sersic pro- 
files with n = \ rather than n = 4 so that the above relation 
between a and 7 does not apply anymore. We have there- 
fore recomputed the mass integral over the same grid, but 
setting now n = 1. This gives us : 

(o, 6, c, d, e, /) ~ (-1.641, 0.644, -0.927, -0.463, -2.491, 0.662) 
thus giving : 

7 ~ (0.609 - 0.736 X a)"^ . 

Unfortunately, setting 7 = 0.14 and solving for a. we get 
a ~ —8.9 so that we are unable to escape the same problems 
as with n = 4. Moreover, our approximating expression (30) 
now works well only over a limited range in {Ceff,ria.) with 
significantly worse residuals ((AX/X) ~ 6%, {/SX/T)rms — 
14%). We therefore caution the reader to not overrate the 
validity of the 7- a relation for n = 1 Sersic models. 



6 A DARK SIZE -MASS RELATION 

Although our main aim was to investigate the size - mass re- 
lation for the visible component of ETGs, it is worth stress- 
ing that the above procedure is quite general. As an inter- 
esting application, we consider dark haloes since it has been 
proposed (Merritt et al. 2006) that the PS model best ap- 
proximates the density profile of dark haloes coming out 
from dark matter only simulations. Actually, the PS param- 
eters change from one case to another so that a general rule 
could not be extracted. Nevertheless, the value of n is not 
too dispersed so that we will repeat the above calculation 
setting n to the average value n = 3.58. As a result, we 
obtain for the slope 7 of the R^f f-M relation : 

7 ~ 1.915 - 0.263 X a . 

An important remark is in order here. Since we are now using 
the PS model for the dark halo, Reff and M can no more 
be estimated unless one fits the model to kinematical (such 
as the velocity dispersion profile) or Icnsing (e.g., Einstein 
rings) data. As such, we cannot provide any estimate neither 
for a or 7. However, it is convenient to reparameterize the 
size - mass relation for haloes in terms of the concentration 
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parameter and the virial mass. Without loss of precision, we 
can simply identify the halo virial mass with the PS total 
mass. In order to define a concentration for the PS halo, we 
first introduce the scale radius Vs defined by the condition : 

dlnp 

dim 

in close analogy with the meaning of in the NFW model. 
Using Eq.(3), we easily get : 



= -2 



n{2 - pn 

bn 



R. 



eff 



(37) 



The concentration cps may be then defined as cps = 
Rvir/vs, and, using R^ir oc mI{^ and i?e// oc MJ.^, we 
finally get : 

CPS oc M-J°" 

with : 

70M = 7 — 1/3 '. 



1.085 + 0.263 X a 



3(1.915 - 0.263 X a) ' 

For < a < 1, wc got 0.19 < "foM < 0.27, so that the con- 
centration is a mild decreasing function of the virial mass. 
Such a behavior is in agreement with the results of Simula^ 
tions, but a direct comparison is not possible because of the 
use of different halo models. 



7 CONCLUSIONS 

In the present paper, we studied the effect of injecting an- 
gular momentum in a Sersic isotropic profile. Our starting 
motivation was to investigate the role of angular momen- 
tum in explaining the discrepancy between FP projections 
(in our case, the Ref / - M* relation) and observations. More- 
over, such a work can also shed light on the differences be- 
tween dwarf and normal ellipticals as evidenced by the same 
FP projections. 

Injecting angular momentum makes the velocity disper- 
sion tensor anisotropic so that the DF has to be accordingly 
modified. The Osipkov-Merritt parameterization provides a 
valuable tool to infer the corresponding DF and hence the 
modified density profile. As a consequence, the size - mass re- 
lation turns out to be changed with respect to the isotropic 
case thus tilting the log Ref / - log M* relation with respect 
to the isotropic value. As we have shown, the slope 7 of 
the size -mass relation is an easy function of the anisotropy 
parameter a. For massive (M* > 2 x 10^° MqETGs, 
our model may be reconciled with the data provided the 
anisotropy profile is tangential in the inner regions (i.e., for 
r/ra< 0.6) to become then radial elsewhere. It is worth not- 
ing that, while a is the only parameter controlling the slope 
of the size- mass relation, all the three quantities (a, Va, Lo) 
enter in determining the zeropoint of the correlation. As a 
consequence, while a can be determined by 7, we cannot 
infer a unique value for the anisotropy scale Va/Reff unless 
we estimate somewhat the reference angular momentum Lq. 

Actually, the unique way to break this degeneracy re- 
lies on fitting the velocity dispersion profile for a statistically 
meaningful sample of ETGs in order to determine, on a case - 
by- case basis, the value of {a,ra/Reff)- Such a test is also 
mandatory in order to check the validity of the assumed 
OM parameterization of the anisotropy profile. Moreover, 



the distribution of the {a,ra/Reff) values can provide an 
estimate of the expected scatter in the log i?e// - log M» re- 
lation thus giving a further cross check of the model. 

It is worth wondering whether our results can help in 
elucidating why dwarfs and giant ellipticals present surpris- 
ing structural differences. Indeed, while dEs appear as the 
natural low luminosity counterpart of giant Es, they do not 
follow the same scaling relations as the latter ones. For in- 
stance, the Ref f - M* relation is rather shallow with a slope 
7 = 0.14 significantly smaller than the one (7 = 0.56) of the 
high mass ellipticals. To this end, we first remember that, 
while both normal and dEs formed from the gravitational 
collapse of primordial density fluctuations, their structural 
and stellar population properties depend on the ability of 
the baryons in a given overdense region to cool and form 
stars. In particular, a feedback mechanism that transfers 
energy back to the interstellar medium is usually invoked to 
explain the dEs structure. For instance, models that invoke 
the cessation of star formation by supernova driven winds 
provide a plausible explanation for the variation of density 
(surface brightness) and metallicity (color) with luminosity 
(Larson 1974; Saito 1979; Vader 1986; Dekel & Silk 1986; 
Arimoto &: Yoshii 1987). On the other hand, in Del Popolo 
(2002) and Del Popolo et al. (2005), it has been shown that a 
similar role may also be played by the acquisition of angular 
momentum. Indeed, the acquisition of a higher amount of 
angular momentum by a less dense region helps contrasting 
the gravitational collapse and hence prevents the cooling of 
baryons and the star formation. 

Based on this latter model, we therefore expect thart 
injecting angular momentum in the dEs DF produces a de- 
crease of mass infall towards the centre thus giving rise to a 
structure with lower central surface brightness. In order to 
qualitatively confirm this picture, it is worth remembering 
that dEs present roughly three different morphologies with 
most bright dwarfs {Mb < ^16) presenting a distinct lu- 
minosity spike in their centre (referred to as the nucleus). 
Lower luminosities dwarfs do not present such a nuclear re- 
gion and may be divided in SO -like and elliptical - like with 
this latter tipycally having a very low surface brightness. 
The different morphology with luminosity can be explained, 
in the previous scheme, as follows: bright dwarfs are born 
from higher density peaks, less subject to tidal torque, and 
as a consequence one expects that more mass freely fall to 
the centre giving rise to the central nucleus. Low surface 
brightness types are born from lower density peaks which 
suffer a larger tidal torque. This produce galaxies without 
a central nucleus, more extended and having low surface 
brightness. Wc have previously seen that in order to ex- 
plain the Reff - Af* scaling law for dEs we need a negative 
value of a implying a large angular momentum term in the 
anisotropic DF, somehow in agreement with this qualitative 
picture. Moreover, the finding that very different values for 
Q are needed in order to explain the reproduced the ob- 
served size -mass relation for dEs and normal Es confirms 
our notion that they have a very different dynamical origin. 

While the outcome of this simplified investigation goes 
in the right direction, more work is needed in order to put 
these preliminary results on a firmer ground. As a first major 
improvement, one has to give off the hypothesis of spherical 
symmetry allowing for rotational flattening. This is partic- 
ularly important since some evidences argue in favour of 
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dEs being the final evolutionary state of tidal dwarfs galax- 
ies (Gentile et al. 2007) thus originating from rotationally 
supported gas rich tidal arms. It is therefore worth investi- 
gating whether injecting angular momentum in the DF of a 
flattened PS model may allow us to still preserve our pre- 
dicted Ref / - Mi, relation and, at the same time, solving the 
problem of the negative a needed to reproduce the observed 
slope. Should this be the case, we could have a strong ev- 
idence in favour of angular momentum playing the leading 
role in determining the structure of dwarf ellipticals. 
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